Number of loops of size h in growing scale-free networks. 
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The hierarchical structure of scale-free networks has been investigated focusing on the scaling 
of the number Nh{t) of loops of size /i as a function of the system size. In particular we have 
found the analytic expression for the scaling of Nh (t) in the Barabasi- Albert scale-free network. We 
have performed numerical simulations on the scaling law for Nh{t) in the BA network and in other 
growing scale free networks, such as the bosonic network and the aging nodes network. We show 
that in the bosonic network and in the aging node network the phase transitions in the topology 
of the network are accompained by a change in the scaling of the number of loops with the system 
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Natural and social systems ranging from protein- 
protein interactions [1] to the World Wide Web [2] can 
be represented as scale- free (SF) networks [3-6], that is, 
as a set of nodes connected by links with special statisti- 
cal properties. These systems are characterized by large 
fluctuations in the individual degrees, i.e. in the number 
of links pointing or leaving each node. In a SF network, 
the statistical distribution P{k) of the degree fc is a power 
law P{k) ^ . Accordingly, if 7 < 3, the individual 
degree distribution displays infinite variance in the limit 
of infinite network size. 

Besides this striking property, models traditionally de- 
fined on regular lattices display peculiar effects when de- 
fined on a scale-free topology instead [7-12]. Such phe- 
nomena have increased the interest in finding universal 
mechanisms to generate SF networks, which seem to arise 
independently in so diverse contexts. 

Recently, these properties have been reproduced by the 
Barabasi- Albert (BA) dynamical model of a random net- 
work. Its simple and natural algorithm, based on a "pref- 
erential attachment" rule [2, 13], triggered an avalanche 
of research activity in the field and generated a rich "zo- 
ology" of SF network models sharing the same fundamen- 
tal ingredients. Though recently alternative mechanisms 
[14,15] for the generation of scale-free network have been 
proposed, in this work we restrict our attention to grow- 
ing scale-free networks with "preferential attachment" . 
But the degree distribution is not the only relevant topo- 
logical quantity. For the characterization of networks it 
is necessary to look at the motifs [16-18] recurrent on it 
and in particular at the number of loops of length larger 
then three [19]. In this work we derive a general formula 
for the number of loops of size h that generalizes the 
known result concerning triangular loops [20] (ft, = 3) for 
a BA network. In particular, we find that the number of 
loops Nh [t) of size h at time t in a B A network scales as 
Nh{t) ^ ("^/2 log(t))'^. In addition, we show that our for- 
mula is robust when tested on other growing SF network 
models, such as the bosonic network model (BN) [21] and 
the aging nodes model (AN) [22]. In the bosonic net- 
work the scaling of the number of loops with the system 
size changes below the Bose- Einstein phase transition be- 



coming a power-law scaling of the type Nh{t) ^ t^. 

The BA model [2,13] was the first and simplest algo- 
rithm generating scale-free undirected networks. In this 
model, a new node is added to the network at each time 
step, and it is connected by a fixed number of links m 
to highly connected existing nodes (preferential attach- 
ment). According to this rule, the probability 11; that an 
existing node i at time t acquires a new link is assumed 
to be proportional to its degree ki(t). This reads 
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The model can be analyzed by a mean field approxima- 
tion [2]. By this approach, one finds that the average 
degree of a node i that entered the network at time ij 
increases with time as a power-law 
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A network built in this way displays a power law degree 
distribution P(/c) ~ k^^. In order to investigate the in- 
homogeneous topology of the network, we define a loop 
of size h (a /i-loop) as a closed path of h links that visits 
each intermediate node only once. 

In the case m > 2, the BA scale-free network is a very 
compact network, with loops of any size. As the network 
evolves, new loops are introduced in the network. By def- 
inition, new loops include the newly added node: indeed, 
a new /i-loop is formed if the new node is connected to 
two nodes already connected by a self-avoiding path of 
size ft, — 2. We indicate with pi^k the probability that the 
nodes i, k, attached to the network at time ti, tk, are con- 
nected by a link. The rate at which new loops of length 
ft are formed is given by the probability that the new 
node k is linked to two existing nodes i and j times the 
probability Pl^~'^it) that they are already connected by 
a self-avoiding path of size ft — 2. Therefore, we write the 
following rate equation for the average number of ft-loops 
Nhit) 
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where the factor i takes into account that each pair of 
nodes i, j has been counted twice in the sums. 

On the other hand, two nodes belong to a /i-loop if 
there is a hnk between them and, besides it, they are 
connected by a self-avoiding path of length h — 1. Let 



dt 



1 /m\3 



(?) 



1 1 1 

'TiTjt 



P^r\t) be the probability that this path exists; 'thus, Integrating (10) we find, in agreement with [20], 



the probability that the link between node i and j be- 
longs to a /i-loop is given by PijPiJ^{t). We obtain the 
average number < Nh{t) > of /i- loops in the system times 
2h by summing this quantity over all the nodes i, j in the 
network. In fact, each loop has been counted 2h times, 
because there are h nodes in the loop, and two possible 
directions. Hence, we can write 
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By replacing Eq. (2) in Eq. (1), the probability that the 
node i is attached to node k is 



Pi.k = 
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Moreover, the probability that a node k is connected with 
nodes i and j is proportional to the probability that i and 
j are already connected, i.e. 



1 

Pi,kPj,k = '^'j^Pi^j- 



(6) 



By replacing this result in (3) and by the definition (4), 
we obtain 

= if- 'X > ^ p, 

Consequently, the rate at which new loops of size h are in- 
troduced in the system is proportional to the mean num- 
ber of loops of size h — 1. 

Equation (7) has a recursive structure that allows its 
integration without any detailed information about the 
probabilities Pl^j{t). In fact, the rate at which new loops 
of length h are formed can be expressed only in terms of 
the number of loops of minimal size (i.e. h = 3), 
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with ( = ^log(i). The number of triangular loops 
< -^3(0 > can be computed directly, for the triangu- 
lar loops increase in time following (3), 
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Since Pij is given by eq. (5), approximating sums by 
integrals we write the rate equation in the form 
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FIG. 1. Scaling of N^it) as a function in a BA network for 
t wp to t = 10^. The Inset shows the values of the measured 
exponent ip defined in(16) for h = 3,4, 5. 

Using eq. (11) in eq. (8), we compute the number of 
loops of size h, < Nh{t) > and we find 

<NH{t)> = [j\og{t)\\l + 0{C')). 
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The expression for the scaling of Nh{t) with the sys- 
tem size t in a BA network does not suggest a practi- 
cal way to measure Nh{t). To this purpose, one has to 
study the symmetrical adjacency matrix a of the net- 
work, whose generic element aij is defined by Uij = 1 if i 
and j are connected and aij = otherwise. Knowing this 
matrix, one directly measures the number of paths start- 
ing from a node i and returning on it after h steps that 
visit intermediate nodes only once. According to this ar- 
gument, the term Nh,{t) has a dominating term of the 
type Si(o'')i,i/(2^) and subdominant terms excluding 
all trivial contributions coming from paths intersecting 
on themselves. Let us assume that the network does not 
contain self-loops, i.e. an = for all i of the network. In 
this case, for /i = 3 we simply have 
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For h = 4 and h 
to show that 



5, by simple arguments it is possible 
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Using relations (13), (14), (15), we can directly measure 
Nh{t) for /i = 3, 4, 5 in the BA scale-free network model 
to check our analytical results. 

In Fig. 1 we show that the scaling of Nh{t) for a BA 
network with m = 2 follows 
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The measured effective exponent ip{h) reported in the 
Inset of Fig. 1 grows with h as expected. Neverthe- 
less, ip{h) differs from the predicted asymptotic behavior 
ip{h) = h in a significant way. This can be explained by 
the fact that we are considering networks of size up to 
t = 10* nodes and we are still far form the asymptotic 
behavior i — > oo. 

To check the robustness of the scaling relations of 
Nfi{t), we have measured the number of loops of size 
h = 3,4,5 in two alternative growing SF network mod- 
els: the bosonic network (BN) [21] and the aging nodes 
network (AN) [22]. 

In the BN bosonic network, each node i is assigned 
an innate quality, represented by a random 'energy' ej 
drawn form the probability distribution p{€i). The at- 
tractiveness of each node i is then determined jointly by 
its connectivity fcj and its energy e^. In particular, the 
probability that node i acquires a link at time t is given 

by 



(17) 



i.e. low energy, high degree nodes arc more likely to ac- 
quire new links. The parameter /3 = 1/T in llj tunes 
the relevance of the quality with respect to the degree 
in the acquisition probability of new links. Indeed, for 
T ^ 00 the probability Ilj does not depend any more 
on the energy and the BN model reduces to the BA 
model. On the other hand, in the limit T ^ only the 
lowest energy node has non zero probability to acquire 
new links. Moreover, in Ref. [21] it has been shown that 
the connectivity distribution in this network model can 
be mapped on the occupation number in a Bose gas. Ac- 
cording to this analogy, one would expect a corresponding 
phase transition in the topology of the network at some 
temperature value Tc- 

In fact, for energy distributions such that (p(e) 
for e — > 0); one observes a critical temperature Tc- For 
T > Tc the system is in the "fit-get-rich" (FGR) phase, 
where nodes with lower energy acquire links at a higher 



rate that higher energy nodes, while for T <Tc a. "Bose- 
Einstein condensate" (BE) or "winner-takes-all" phase 
emerges, where a single nodes grabs a finite fraction of 
all the links. We simulated this model assuming 



p{e) = {e + l)e^ and e e (0, 1) 
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where 9 = 1. For this distribution, it has been shown 
that Tc{9 = 1) = 0.7 [21]. In particular, in Fig. 2 we 
report Nh{t) for /i = 3,4, 5 as a function of the number 
of nodes t in the network, at temperature T = 1.5 and 
T = 0.5 respectively. 

According to simulations, in correspondence with the 
phase transition, the scaling of Nfi{t) drastically switches 
from a BA-like behavior 
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to a power-law scaling 

Nh{t) (X f^e*) for T < T^. 
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FIG. 2. Scaling of the number of loops Nh{t) of size h with 
the system size t in the FGR and in the BE phase of a bosonic 
network. The plotted data have been obtained for a network 
with p(e) = 2e and e £ (0,1) (which has Tc = 0.7 [19]) at 
T = 1.5 (FGR piiaso) and at T = 0.5 (BE phase). In the 
Insets we report the value of the exponents ij) (2.6,3.9,5.2) 
and ^ (0.78,1.28,156) found in the simulations for /i = 3, 4, 5 
respectively. 

We have measured the scaling of Nh{t) with the sys- 
tem size t also for a growing network with aging nodes 
(AN) introduced in [22] . The model has been motivated 
by the observation that in many real networks, e.g. the 
scientific citations network, old nodes are less cited than 
recent ones. In the goal of representing this feature, the 
probability 11^ to attach a link to a node i arrived in the 
network at time U is modified to be 
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where a is an external parameter. As in the BA model, 
a new node is connected to m existing nodes. The re- 
sulting structure of such a network strongly depends on 
the constant a. For a < 1, the degree distribution is a 
power law P{k) ~ k~'^ with an exponent monotonically 
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increasing from 7 = 2 in the limit a— oot0 7^ooin 
the hmit a 1: on the other hand, for a > 1 no power 
law is observed in the degree distribution. Therefore, this 
model reproduces a SF network only in the region a > 1. 

Moreover, as observed in [22], in the limit a —00 
the oldest node is connected to an increasing fraction of 
all the links, reminding the condensation observed in a 
bosonic network. To take into account this phenomenon, 
we introduce a value a* such that for a < a* the fraction 
of links attached to the most connected node exceeds a 
finite threshold F. We expect the scaling of Ni^{t) to 
be different in the three regions a > 1, a G {oi*,l) and 
a < a*. 



of the logarithm of the system size. Moreover, we have 
observed that indeed this scaling seems to be a marking 
feature of growing scale-free networks with preferential 
attachment. In particular, topological phase transitions 
in the bosonic network and in network including aging 
nodes are accompanied by a drastic change in the scaling 
of Nh{t) with the system size. 

The authors arc grateful to Yi-Cheng Zhang and 
Alessandro Vespignani for their useful comments. This 
work has been partially supported by the European Com- 
mission - Fet Open project COSIN IST-2001-33555. 
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FIG. 3. The order parameter kmax/mt for a network with 
aging of the nodes, size t = lO** and m = 2. We distinguish 
between three region of the phase spaceia > 1, a £ (—1, 1) 
and a < — 1. In the Insets we report the typical behavior of 
Nh{t) as a function of t for h = 3, 4, 5 in the three regions. 

We measured the total fraction of links kmax/n^t at- 
tached to the oldest node in a network with m = 2 and 
t = 10^ nodes. The threshold has been fixed at = 0.1, 
in order to distinguish the 'condensate' phase from the 
simple scale-free phase. The value for a* was found to 
be a* = —1. We then measured the number of loops of 
size /i = 3, 4, 5 for networks made of up to 10^ nodes in 
the three ranges of value of a. We have observed that, 
for a > I, the number of loops of size h scales linearly 
with t (at least for /i = 3, 4, 5). In the Inset (a) of Fig. 3 
we report the data for a = 1.5. On the contrary, for 
a € {—a*, 1) we measured the scaling 

ylog(i)) 

with ip{h) a monotonic function of h. 
Fig. 3, data for a = 0.5 are reported, 
region a < a* = — 1, Nh{t) becomes proportional to a 
power-law of the system size, as it is shown in Inset (c) 
of Fig. 3, referring to the case a = —5. 

In conclusion, we have introduced the number of h- 
loops Nh{t) for a network of t nodes as a characterizing 
quantity for random networks. We have observed that in 
the BA scale- free networks Nh,{t) scales as a power law 
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Finally, in the 
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